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We study the spectrum of low-lying eigenmodes of the kinetic operator for scalar particles, in the color adjoint 
representation of Yang-Mills theory. The kinetic operator is the covariant Laplacian, plus a constant which 
serves to renormalize mass. In the pure gauge theory, our data indicates that the interval between the lowest 
eigenvalue and the mobility edge tends to infinity in the continuum limit. On these grounds, it is suggested that 
the perturbative expression for the scalar propagator may be misleading even at distance scales that are small 
compared to the confinement scale. We also measure the density of low-lying eigenmodes, and find a possible 
connection to multi-critical matrix models of order m — 1 . 



I. INTRODUCTION 

In recent years it has been recognized that kinetic operators 
in confining gauge theories may have a low-lying spectrum of 
localized eigenstates [1-3]. By "kinetic operator" we mean, 
e.g., a Euclidean Dirac operator or covariant Laplacian oper- 
ator in a background gauge field, with the possible addition 
of a constant representing a mass term. The "mobility edge" 
is a point in the spectrum between an interval of spatially lo- 
calized eigenmodes, and the bulk of eigenmodes which are 
extended over the full volume. 1 

In condensed matter physics, it has long been known that 
the Hamiltonian of an electron moving in a stochastic poten- 
tial has a low-lying spectrum of localized eigenstates [4]. If 
the energy at the mobility edge is higher than that of the Fermi 
surface, then the material is an insulator [5]. This is because 
the propagation of a wavepacket is essentially an interference 
effect among energy eigenstates which are spatially extended. 
If there are no available extended modes, then there is no par- 
ticle propagation through the medium. 

This condensed matter example motivates us to ask the fol- 
lowing question: In a confining lattice gauge theory, what 
is the interval in GeV (Dirac operator) or (GeV) 2 (covariant 
Laplacian) between the lowest eigenvalue of kinetic operator, 
and the mobility edge, as we take the continuum limit? If this 
interval shrinks to zero in physical units, or is populated by 
only a finite number of states in the continuum limit, then the 
relevance of localized states to particle propagation is ques- 
tionable. On the other hand, suppose the energy interval of 
localized states goes to infinity, in the continuum limit, for a 
kinetic operator of some type. This fact may then have rad- 
ical implications for gauge theories with matter fields of the 
corresponding type. 

In this article we investigate the low-lying spectrum of the 
covariant Laplacian in the adjoint representation of the SU(2) 



gauge group. Up to a constant, the covariant Laplacian is the 
kinetic operator for adjoint scalar fields. The adjoint repre- 
sentation is chosen because we know, from the work in ref. 
[2], that there is already something odd about the spectrum: 
the lowest eigenmode is localized in a region whose volume 
goes to zero in physical units (but infinity in lattice units) in 
the continuum limit. The consequences of this fact for the 
propagation of scalar particles was not so clear. In our present 
paper we extend our study to the full spectrum of localized 
states. Our objective is to (i) find how the average interval 
AX mo b between the lowest eigenvalue and the mobility edge 
varies with coupling; and (ii) compute the density of eigen- 
modes in this interval. With those results in hand, we discuss 
possible consequences for the two-point function of adjoint 
scalar fields in confining gauge theories. 



II. SCALING OF THE LOCALIZATION INTERVAL 

The kinetic operator for scalar fields, in the adjoint repre- 
sentation of the SU(2) gauge group, is given by 
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is the covariant lattice Laplacian, with link variables 
Uf = \Tr[o a U,{x)o b Ul(x)] 
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1 In fact there is usually an interval of localized states at both the lower and 
the upper ends of the spectrum, and therefore two mobility edges. 



in the adjoint representation, and ;«q may be taken negative. 
We compute numerically the low-lying eigenvalues and eigen- 
modes of the covariant Laplacian operator 

-flgtfCv) (2-4) 

in thermalized lattice configurations, and from these we com- 
pute the inverse participation ratio (IPR) of each eigenmode 

/ra„ = v£lCto<MI 2 (2.5) 
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From our previous work [2] we know that there is some range 
of eigenvalues of the covariant Laplacian, lying between the 
lowest eigenvalue ko and the mobility edge k mo b, whose cor- 
responding eigenmodes are localized, and whose IPR's grow 
linearly with lattice volume V. For k n > X mo b the eigenmodes 
are extended, and the IPR's are all of (9(1) at large volumes. 
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FIG. 1: (IPR„)/{IPR ) vs. Ak n = (k n -Xa)/a 2 at p = 2.35, on lat- 
tice volumes from 14 4 to 24 4 . The data points at each volume are 
connected by lines. 
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FIG. 2: Same as Fig. 1, except lines connecting data points at dif- 
ferent volumes are omitted. The straight line is a best fit through the 
combined all- volumes data in the interval AA, G [0 : 0.045] . 

Let a denote the lattice spacing. We are interested in deter- 
mining how the interval of eigenvalues in physical units 

(Knob — k()) 
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varies as coupling p increases, and a — > 0. Our method is 
based on the following observation: If we plot 
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FIG. 3: {lPR n )/{IPR Q ) vs. AX„ = <A„ - A )/a 2 at P = 2.5, on lat- 
tice volumes from 16 4 to 24 4 . The data points at each volume are 
connected by lines. 
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FIG. 4: Same as Fig. 3, except lines connecting data points at dif- 
ferent volumes are omitted. The straight line is a best fit through the 
combined all- volumes data in the interval AA, € [0 : 0.06]. 



at a given coupling p, then the data points at different volumes 
tend to cluster around the same straight line at small AA,, as 
seen in Figs. 1-2, at p = 2.35, and Figs. 3-4 at P = 2.5. There 
does not appear to be a strong volume dependence in the data 
at AA, < Akmot,. Above the mobility edge we would expect 
R„ —> in the infinite volume limit, since in that limit the IPR's 
of localized and extended states tend to infinity and to finite 
values, respectively. Fitting a straight line 



/(AA,) = 1 



Ak 

Akint 
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to the combined volume data at small AA,, the intercept AA. mf 
with the x-axis is a reasonable measure of the width of the 
peak in the R vs. AA, data. Even if AA,,,,; is not precisely the 
same as Ak mo b, we expect these quantities to scale with p in 
the same way. The result, for p e [2.2,2.5] is shown in Table 
I, and we see that the data is consistent with 



Ak in 



0.045(3)GeV 



(2.10) 
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TABLE I: The intercept AA.,„ ( vs. (3, as determined from the best 
linear fit (2.9) to data in the intervals and volumes shown. AXi„ t and 
AXi nt a are in units of GeV 2 and GeV, respectively. 
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2.20 
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[0,0.030] 


0.038(1) 
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2.30 


12 4 ,14 4 ,16 4 ,18 4 


[0,0.040] 
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0.044 


2.35 


14 4 ,16 4 ,18 4 ,20 4 ,24 4 


[0,0.045] 
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0.047 
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14 4 ,16 4 ,18 4 ,20 4 ,24 4 
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0.047 


2.50 


16 4 ,18 4 ,20 4 ,24 4 


[0,0.060] 


0.099(4) 
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Since AXi„ t ~ AX mo b, the implication is that 

limAA. m „ fc =°o (2.11) 

We reserve discussion of this result to section IV, below. 



III. DENSITY OF LOCALIZED STATES 

The covariant Laplacian operator in adjoint representation 
can be thought of as a 3V x 3V random matrix, where the 
factor of 3 comes from the color index. We are interested in 
computing the density of localized eigenstates at the low end 
of the spectrum. Suppose this has the form 

p(X)=K(X-X () ) a (3.1) 

where we normalize the density of states such that 

/ dXp(X)=3 (3.2) 
JXo 

Then the number of eigenmodes with eigenvalues less than X 
is 

n(X)=V / dX p(X) (3.3) 
JXo 

A standard method for determining a is based on the observa- 
tion that under the eigenvalue rescaling 

z = V 1/(1+a) (X-X ) (3.4) 

the corresponding eigenvalue number n(z) becomes indepen- 
dent of V . This is readily verified: the number of eigenvalues 
dn(X) in an interval 8 A. around X is 

8n(X) = Vp(X)8X (3.5) 

Rescale z = V P (X — A-o). Then, using the assumed form (3.1) 
for p(A-), we have 

5n(z)=KV 1 -P( 1+< V5z (3.6) 

Choosing p = 1/(1 +a), 8n(z) is volume independent. So in 
order to compute the density of states, we compute n(X) at a 
variety of lattice volumes, and look for a constant a such that, 



under the rescaling (3.4), the curves for n(z) computed at each 
lattice volume coincide. 

In numerical simulations, we evaluate the first N ev eigenval- 
ues of the covariant Laplacian in a set of N am f independent 
thermalized configurations. The resulting N am f x N ev eigen- 
values of all N con f configurations are then sorted from lowest 
to highest regardless of configuration or eigenmode number. 
Then, if X m is the m-th eigenvalue in the ordered list of eigen- 
values, we identify 

m 

<X m ) = (m = l,2,...,N conf x N ev ) (3.7) 

Nconf 

The maximum value is n(X max ) = N ev , where X max is the 
largest eigenvalue in the list. Other values of n(X) for X < X max 
and X ^ X m are obtained by interpolation (data points are con- 
nected by lines in the figures shown). The numerical results 
for p = 2.5 on L 4 lattices at L = 14, 16,20,24 are shown in 
Fig. 5. 
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FIG. 5: Cumulative number of eigenmodes n(X) with eigenvalue less 
than A, at (3 = 2.5. Data shown is for lattice volumes from 14 4 to 24 4 . 

In Fig. 5 the n(X) curve is different for each lattice volume. 
This is, of course, not surprising, since as volume increases, 
so does the number of eigenvalues in any given interval AA. 
Since we are interested in checking universality under rescal- 
ing, right at the very end of the spectrum, we concentrate on 
the region where n(X) < 10. In Figs. 6-8 we show «(z), un- 
der the rescaling (3.4), for a = 1,2,4 respectively. The data 
seems compatible with universality at a = 2, in which case we 
would have 

n{z) = |z 3 (3.8) 

Fig. 9 shows a best cubic fit (solid line) to the a = 2 scaled 
data at all volumes, in the interval z £ [0,2]; K is determined 
to be 3.8. The cubic curve appears to fit the data at p = 2.5 
quite well, strengthening the case that a = 2 is the correct 
exponent. 

The exponent a = 2 may be of some significance. Accord- 
ing to ref. [6], even-integer exponents a = 2m are obtained 
in simple large-N hermitian matrix models with polynomial 
potentials, when the coupling constants in the potential are 
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FIG. 6: Cumulative eigenvalue number n(z) vs. the scaling variable 
z = yi/(i+«)(X-Xo)atp = 2.5anda=l. 



FIG. 8: Same as Fig. 6, with a = 4. 
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FIG. 7: Same as Fig. 6, with a = 2. 



FIG. 9: Best fit of n(z) = ^Kz 3 to the p = 2.5 data (combined vol- 
umes) at a = 2. 



fine-tuned to achieve criticality. It appears that the low-lying 
spectrum of the covariant Laplacian corresponds to that of a 
multi-critical matrix model of order m = 1 . 

It is also of interest to study the eigenvalue spacing distri- 
bution of the "unfolded" spectrum [7], which has previously 
been studied, in connection with random matrix theory, for 
the Dirac spectrum of SU(2) and SU(3) pure gauge theory. 
The unfolded spectrum is a rescaling of the spectrum such 
that the average spacing between neighboring eigenvalues is 
unity, and a quantity of interest is the distribution of fluctua- 
tions around that average value. The procedure [8] is as fol- 
lows: Let X" be the z'-th eigenvalue of the covariant Laplacian 
in the n-th lattice configuration, with X" > X"_ l . The set of all 
{A,"} of all configurations is then sorted in ascending order, 
and we denote by N" the location (from 1 to N con f x N ev ) of 
X" in the sorted list. For the n-th configuration, the eigenvalue 
spacing s" between neighboring eigenvalues in the unfolded 
spectrum is defined as 



' N f 



(3.9) 



(i.e. regardless of i, n). In the case of the Dirac operator [8], 
the distribution P(s) is well described by a certain Wigner dis- 
tribution. For the spectrum of the covariant Laplacian, our 
numerical results show that eigenvalue spacing distribution of 
the unfolded spectrum is also well described by one of the 
Wigner distributions, namely, the orthogonal distribution 



(3.10) 



Let P(s) denote the probability distribution of these spacings 



as seen in Fig. 10 for P = 2.35 on an 18 4 lattice volume. We 
have also checked that eq. (3.10) gives a good fit to the eigen- 
value spacing distribution at p = 2.4 and p = 2.3. 

Returning to the density of states p (X) of the original spec- 
trum, we find that the quadratic (a = 2) power behavior con- 
tinues well past Xj nt , but does not persist thoughout the spec- 
trum. Fig. 9 showed an excellent cubic fit to the cumula- 
tive n{z) data, but this was for low-lying eigenmodes with 
n(z) < 12. In Figs. 11 and 12 we show all of the data for 
n(z), obtained at p = 2.5 on 16 4 and 20 4 lattices, together 
with the cubic fit to the data at all volumes. The vertical line 
in each graph indicates the value of z which corresponds, at 
each volume, to X — Xo = 0.05. In both cases the data for n{z) 
begins to diverge away from the cubic fit at about this value of 
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FIG. 10: Unfolded eigenvalue spacing distribution for (3 = 2.35, L = 
18. 



A — Xq. In physical units, Ak = (A — Xo)/a 2 at A, — Ao = 0.05 
is roughly two and a half times our estimate for AA,„ r at this 
coupling. 
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FIG. 11: Cumulative eigenvalue number n(z) vs. z at P = 2.5, a = 
2 on a 16 4 lattice volume. The data shown includes all calculated 
eigenmodes, which extend beyond the mobility edge. For reference, 
the vertical line intersects the x-axis at a value of z corresponding to 
X -Xq = 0.05 (see text). 



Given the density of states we can compute the phase-space 
volume, in physical units, taken up by states below the mo- 
bility edge. For purposes of comparison we begin with the 
free case, where it is easy to show that the density of states in 
D = 4 dimensions is p(X) = X/% 2 . The total number of states 
in the eigenvalue interval from X = to X st0 p = Aa 2 is 



X<Aa 2 

X 



= v f 

JO 



dX 



(3.11) 



where a is the lattice spacing, and A is some fixed number 
in units of (GeV) 2 . Now divide by the physical volume, to 
get the number of states per unit volume, T|, which lie in the 
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FIG. 12: Same as Fig. 1 1, on a 20 4 lattice volume, 
eigenvalue interval X/a 2 < A: 
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= ~A f dX ^2 
fl 4 JO It 

1 f A 



i 1a 2 

2 7T 2 



(3.12) 



Then if A is a finite cutoff in physical units, r| is also finite, 
even as a — > 0. 

We make the same computation for localized states of the 
covariant Laplacian, which lie in the eigenvalue interval X £ 

[Xo,X mo h] 



Vphys 

1 ('Knob r. 

= - / dXK(a)(X-X f 
cr Jx Q 

K(fl) (X mob - A-p) 3 

3 a 4 

Denoting eigenvalues in physical units as X phys 
X = X phys a 2 

we have 

K(fl) 2 p hys yphys,3 
T\ - —^~ a y K mob - A ) 

According to the data in the previous section, 

X^-X^=M/a 
where M w 0.045 GeV. This means that 

1 K ( a ) 

ri = — —M 5 
1 a 3 



(3.13) 



(3.14) 



(3.15) 



(3.16) 



(3.17) 
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The question is then how the factor K in the density of states 
(3.1) varies as a function of a. In the fit to the n(z) data at 
P = 2.5, we found K w 3.8. Repeating the same analysis at 
P = 2.35, we find K w 2.8. Thus K does not seem to decrease 
with smaller lattice spacing; if anything, it increases. 

From this data, we draw a remarkable conclusion: The 
number of states below the mobility edge, per unit physical 
volume, diverges in the continuum limit. This means that if 
we compute the contribution to any observable, keeping only 
a finite number of the lowest scalar field eigenmodes per unit 
volume, then all of the contributing states would be localized 
states in the continuum limit. This is an unexpected feature of 
adjoint scalar fields in confining gauge theories, and suggests 
that the particle propagator is dominated, at scales which are 
large compared to the lattice spacing, by the localized eigen- 
modes. In the next section we present a second argument lead- 
ing to the same conclusion. 



IV. LOCALIZATION AND THE SCALAR PROPAGATOR 

In a free theory, the eigenmodes of scalar and fermionic ki- 
netic operators are simply plane wave states, and this is also 
the starting point of the weak-coupling perturbative approach 
to gauge theories. In perturbation theory, it is assumed that the 
eigenmodes of kinetic operators such as the covariant Lapla- 
cian, in a suitable (e.g. Landau) gauge, are approximately 
plane wave states, at least for Euclidean momenta which are 
large compared to Aqcd- On the other hand, we have seen 
in section II that for the adjoint covariant Laplacian, the in- 
terval between the lowest eigenvale A-o and the mobility edge 
X m „b tends to infinity, in physical units, in the continuum limit. 
What this implies is that the contribution of finite momenta 
extended modes to the scalar particle propagator is negligi- 
ble, unless we are willing to tolerate a mass subtraction which 
would introduce tachyonic modes into the theory. 

On the lattice, the scalar particle propagator in the quenched 
(no scalar loop) approximation has the form 



AV, 



phys 



-T&P 4 is the momentum space volume in physical 



G a6 (*-y)=£ 



X n +ml 



(4.1) 



where the VEV is evaluated in the pure gauge theory 
with some appropriate gauge choice. By "finite momen- 
tum extended modes", we mean extended eigenmodes (X > 
'Knob) whose Fourier components are negligible for momenta 
\p\ja < P, where P is a momentum cutoff which is large in 
physical units, but small compared to 1 /a. Consider the con- 
tribution to the scalar propagator, in the continuum limit, due 
to these extended, finite momenta eigenmodes. Denoting this 
contribution by G'(x — y), we have in physical units 



-lab 
7 phys 



X„+ml 



(4.2) 



where Yl n denotes the sum over the finite momenta, extended 
eigenmodes. The number of such eigenmodes cannot ex- 
ceed the number of lattice momenta p satisfying the restric- 
tion \p\/a < P. This number is of order a 4 VAVp hys , where 



units, and V is the lattice volume in lattice units. We can easily 
estimate the magnitude of G' phys from 



£'l ~ a 4 VAvf ys 



so that 



G (x - y) 



AV' p 7hys 



(X' + m 2 )/a 2 



(4.3) 



(4.4) 



where X' is the magnitude of a typical eigenmode in the range 
of eigenmodes summed by Yl n - O n tne other hand, in the con- 
tinum limit, X' /a 2 — ► °°, and therefore the contribution to the 
scalar propagator from finite momentum extended modes is 
negligible, for any > 0. This is by no means special to 
confining gauge theories; it is also true in QED. The resolu- 
tion in QED is to allow for a negative bare mass term < 0, 
adjusted so that A-o + OTq is 0(a 2 ). But for the adjoint scalar in 
Yang-Mills theory, this choice of counterterm is inadequate. 
The problem is that for eigenvalues X„ > X mo \ contributing to 
G'(x-y) 



X n + ml 



> 



Xq - 



■ +AX, 



■mob 



(4.5) 



■n-rm§ = 



Therefore, even if we chose a counterterm such that A-q 
0, the denominator in eq. (4.4) would still be of order AX mo b 
1 1 a, which means that G' —* in the continuum limit. The 
only way to avoid this is to choose a counterterm such that 
Xmob + rriQ is 0(a 2 ). But in that case, the eigenvalues X„ + 
of the kinetic operator for the localized eigenmodes X„ < X mo b 
are negative, i.e. tachyonic. 

In a quenched theory the bare mass constant can be cho- 
sen at will, and tachyon modes in the scalar kinetic term are 
not excluded. However, in a well-defined field theory with 
a dynamical scalar field there can be no true tachyon modes; 
these only appear in perturbative calculations around a false 
vacuum state. Thus the quenched scalar propagator can only 
approximate the scalar propagator in the unquenched theory 
when the lowest eigenmode Xo + of the kinetic operator is, 
on average, greater than zero. 

In D=4 dimensions, SU(2) gauge-Higgs theory with the 
scalar field in the adjoint representation (aka the Georgi- 
Glashow model) is known to have two distinct phases: a con- 
finement phase and a Higgs phase [9] 2 Our quenched calcula- 
tion would be relevant as an approximation to the scalar prop- 
agator of SU(2) gauge-Higgs theory in the confined phase. In 
this theory 



J DU J Dcp < 



-(SyM+Sif) 



2 In contrast, a gauge-Higgs theory with the scalar in the fundamental repre- 
sentation has only one phase, and the asymptotic string tension is vanishing 
for all finite gauge and Higgs couplings (cf. ref. [10]). 
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(4.6) 



where 



S YM [U] 
p -m[u] 



-P£|Tr[f/f/f/f/] 

E [s ( p(- £) ) ( i ) + 3»4 re <p 2 + y<p 4 ] 

Z^t/] = j Dye-WM (4.7) 



The covariant Laplacian and scalar field are in the color ad- 
joint representation; all color indices are implicit. 

Define expectation values in a fixed background gauge field 

as 



(Q) 



and in particular, 



u = i-y"z)q>G[q>,tf]e-*P (4.8) 



f_D£7 ((p"(x)(p fe (>-))t/ e -(fr*+^0 

f £>£/ e -(^M+AS) 



(4.9) 



In the quenched approximation, it is assumed that AS is only 
a small correction to Sym in the confined phase of the theory. 

As before, let (j)„ (x) denote the eigenstates of the covariant 
Laplacian for fixed U, and expand the scalar field 



Defining 



and assuming 



— = (*„*„) t 



(®r&m)u ~ for m^n 



we have 



<<p»<p"M)=E 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



The quenched approximation is based on the assumption 
that it is possible to choose some values for p and in the 
quenched theory, dependent on the gauge coupling and the 
Higgs couplings m 2 harel yin the confined phase of gauge-Higgs 
theory, such that 



A„ 



gauge-Higgs 



' M X„ + m2 /, 



'0 / pure YM 

(4.14) 

where the expectation values on the lhs and rhs are taken in the 
unquenched and pure Yang-Mills theories respectively. The 
argument for the validity of this approximation is the usual 
one: large-N, and the general expectation that (i) as long as 
the system remains in the confined phase, the effect of matter 
loops on the vacuum state is not very drastic; and (ii) the main 
effect of the ycj) 4 term on the scalar particle propagator is to 
renormalize the mass term. 



In gauge-Higgs theory, the A„ must, from their definition, 
be positive semi-definite; the quenched approximation (es- 
sentially (A„) « (k„ +ml)) can only be relevant to the un- 
quenched theory for (k n + OTq) > 0. In that case, as we have 
seen, the contribution of finite momentum extended modes 
to the scalar propagator is negligible, and it is the localized 
modes which dominate scalar particle propagation at distances 
which are large compared to the lattice scale. This would 
mean that ordinary weak-coupling perturbation theory goes 
very wrong for adjoint scalar particles, even at distance scales 
which are quite small compared to the confinement scale. 3 



V. CONCLUSIONS 

In a previous article [2] we had found something odd about 
the spectrum of the covariant Laplacian in the adjoint repre- 
sentation: the lowest eigenmodes appear to be localized in a 
volume which shrinks to zero, in physical units, in the contin- 
uum limit. In the present work we have extended our study 
to the full interval of localized states, and have found other 
surprising features: First, the range of eigenvalues AX mo b of 
the localized eigenmodes tends to infinity, in physical units, 
in the continuum limit. Secondly, the density of eigenmodes 
rises from zero quadratically, up to the mobility edge (sug- 
gestive of a connection to multi-critical matrix models of de- 
gree m—1), and the number of localized eigenmodes per unit 
physical volume is infinite in the continuum limit. We must 
add a caveat that these conclusions rest on the results of nu- 
merical simulations of SU(2) lattice gauge theory, carried out 
at couplings between p = 2.2 and P = 2.5, and lattice volumes 
up to 24 4 . It is, of course, not excluded that the trend in the 
data could change at higher couplings and/or larger volumes. 

An infinite range of eigenvalues between the lowest eigen- 
value and the mobility edge has an interesting consequence. If 
mass counterterms leading to tachyonic modes are excluded, 
then the quenched scalar particle propagator is dominated by 
localized modes; there is a negligible contribution from ex- 
tended eigenmodes corresponding to finite physical momenta. 
Exclusion of tachyonic modes is necessary if the quenched 
propagator is to be a reasonable approximation to the scalar 
propagator in a gauge theory with dynamical scalar fields. 
Again there is a caveat: We do not really know if an infinite 
range of localized eigenmodes persists in the unquenched the- 
ory. This will require further, computationally more intensive, 
investigations. 

If the propagator for adjoint scalar fields is completely dom- 
inated by localized states, when evaluated in the confined 
phase of gauge-Higgs theory in some appropriate gauge, then 
this would raise some doubts about the validity of perturbation 
theory, at least as applied to adjoint scalar fields. Although 
one naturally expects weak-coupling perturbation theory to 



' Another short-distance phenomenon which is missed by perturbation the- 
ory is the excess lattice action found on thin P-vortex sheets and monopole 
lines; c.f. ref. [11]. 
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break down at a distance scale comparable to the confinement 
scale, it is generally believed that this procedure should pro- 
vide correct answers for short-distance quantities. If, however, 
the scalar propagator is dominated, even at short distances, by 
localized eigenmodes, then weak-coupling perturbation the- 
ory may be misleading. This is an interesting (and obviously 
radical) possibility, which calls for further investigation. 
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APPENDIX: ALTERNATIVE FITS 



of X (which is divided into ten intervals) depends on the high- 
est and lowest values of X n found in all lattice configurations. 
This means that lattice configurations in which the lowest 
eigenvalues are well below the average may introduce data 
points in the graph whose lowest X value is also well below 
the average Xq. Nevertheless, a linear fit of IPR l l p vs. X, with 
p^l, might conceivably be superior to a linear fit of I PR vs. 
X in the neighborhood of the mobility edge, and this in turn 
would give a better estimate for X. mo ^. We have concluded, 

0.2 | 1 1 1 1 1 1 1 1 1 1 




In addition to the linear fit of IPR to X, we have explored 
some other possible fits and fitting procedures for estimating 
X mo b. In particular we have tried to fit various fractional pow- 
ers of the IPR to X, and to bin the data in different ways. As 
an example, we have subdivided the interval between minimal 
and maximal value of X into 10 subintervals and calculated the 
average IPR 1 '''' of the eigenvalues in each interval. We asso- 
ciate, with each average IPR 1 /'', a value of X in the center of 
the corresponding X interval, and then fit these data points to 



=AX + B 



IPR 

I? 



The x-axis intercept A.,, 



Xint — — 



B 



(A.l) 



(A.2) 



is taken to be an estimate of the mobility edge X m( ,b ~ Xi„ t . 

A plot of (IPR/L 4 ) 1 / 4 vs. X at P = 2.4, on a 20 4 lattice is 
shown in Fig. 1 3 ; the data obtained at p = 2 .4 on 1 6 4 , 20 4 , 24 4 
lattice volumes are displayed together in Fig. 14. We should 
note that these plots tend to underestimate A-o, in that the range 
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FIG. 14: Dependence of (2£ 



0.25 



on X for various L and P = 2.40. 



however, after some experimentation with different fits, that 
the data is not really adequate to convincingly determine the 
optimal value of p near the mobility edge. 
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